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Abstract 

In  this  study,  we  explain  and  quantify  the  dynamical  causes  of  striking  differences  in  the  acoustic 
transmission  data  collected  on  the  shelf  and  shelfbreak  in  the  north-eastern  Taiwan  region  within  the 
context  of  the  Quantifying,  Predicting  and  Exploiting  (QPE)  uncertainty  2008  Pilot  Experiment.  To 
do  so,  we  employ  our  own  coupled  oceanographic(4D)-acoustic(Nx2D)  modeling  systems  with  ocean 
data  assimilation  and  a  best-fit  depth-dependent  geo-acoustic  model.  Eor  the  first  time,  predictions  are 
compared  to  the  measured  acoustic  data,  showing  skill.  We  also  study  the  sensitivity  of  our  results  to 
uncertainties  in  several  factors,  including  geo-acoustic  parameters,  bottom  layer  thickness,  bathymetry 
and  ocean  conditions.  Our  simulations  first  reveal  that  the  lack  of  signal  received  on  the  shelfbreak  is 
due  to  a  dramatic  20  dB  increase  in  transmission  loss  (TL)  caused  by  bottom  trapping  of  sound  energy 
during  up-slope  transmissions  over  the  complex  and  deeper  bathymetry.  Sensitivity  studies  on  sediment 
properties  show  larger  but  isotropic  TL  variations  on  the  shelf  and  smaller  but  more  anisotropic  TL 
variations  over  the  shelfbreak.  Diverse  thicknesses  of  sediments  lead  to  only  limited  effects  on  the  TL. 

The  small  bathymetric  data  uncertainty  are  modeled  and  also  shown  to  lead  to  small  TL  variations. 

We  discover  that  the  initial  transport  conditions  in  the  Taiwan  Strait  can  affect  acoustic  transmissions 
downstream  more  than  100  km  away,  especially  above  shelfbreak.  Simulations  also  reveal  internal  tides 
and  quantify  their  spatial  and  temporal  effects  on  the  the  ocean  and  acoustic  fields.  One  type  of  waves 
found  are  shelfbreak  internal  tides  propagating  up-slope  with  wavelengths  around  40-to-80  km  and 
vertical  displacement  of  isotherms  of  lO-to-20  m.  These  waves  lead  to  variations  of  broadband  TL 
estimates  over  5-6  km  range  that  are  more  isotropic  and  on  bearing-average  larger  (up  to  5  dB)  on  the 
shelf  than  on  the  complex  shelfbreak  where  the  TL  varies  rapidly  with  bearing  angles. 

Index  Terms 

Data  assimilation,  underwater  acoustics,  shelf  and  shelfbreak  sound  propagation,  uncertainty  quan¬ 
tification,  interdisciplinary  modeling 

I.  Introduction  and  Motivation 

As  one  of  the  major  applieation  of  underwater  aeousties,  sonar  performanee  predietion  requires 
modeling  the  aeoustie  field  evolution.  In  the  littoral  oeean  environment,  the  time  and  spaee  seales 
relevant  for  sueh  predietions  ean  be  minutes  to  a  week  and  hundreds  of  meters  to  tens  of  kilo¬ 
meters.  The  parameters  inelude  the  four-dimensional  oeean  and  seabed  fields.  They  are  eomplex 
to  prediet  and  ean  have  signifieant  uneertainties  [1],  [2].  Methods  and  systems  that  foreeasts 
the  oeean,  the  seabed  and  the  aeousties  in  an  integrated  fashion  have  only  been  developed  and 
utilized  reeently  [3]-[6].  Aeoustie  data  were  previously  not  available  for  evaluations  of  sueh 
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systems.  An  objective  of  the  present  investigation  is  to  compare  predictions  of  such  coupled 
systems  to  in  situ  acoustic  data.  Even  though  important,  such  comparisons  had  not  been  done 
before. 

Our  approach  is  based  on  coupling  data-assimilative  environmental  and  acoustic  propagation 
models  with  ensemble  simulations,  as  developed  by  [4],  [5],  [7],  [8].  This  approach  was  applied 
with  adaptive  sampling  and  onboard  routing  within  the  Focused  Acoustic  Forecasting-05  (FAF05) 
excercise  [9].  During  the  PN07  exercise  in  Dabob  Bay,  acoustic  transmission  loss  fields  were  also 
coupled  to  data-assimilative  ocean  fields  by  Xu  et  al.  [10]  so  as  to  characterize  impacts  of  wind 
forcing  and  tidal  forcing  on  acoustic  fields  and  performance.  Within  the  scope  of  Battlespace 
Preparation  2007  (BP07),  Fam  et  al.  [11]  coupled  acoustic  and  ocean  models  at  sea  in  real¬ 
time,  Rixen  et  al.  [12]  utilized  super-ensemble  prediction  techniques  for  acoustic  inversion  and 
tomography,  Carriere  et  al.  [13]  investigated  full-field  tomography  and  tracking,  and  Martins 
and  Jesus  [14]  studied  acoustic  prediction  as  a  Bayesian  estimation  problem.  The  New  England 
Shelbfreak  front  region  e.g.  [15]-[17]  and  the  shallow  Asian  Seas  e.g.  [18]-[22]  have  also  seen 
major  ocean  and  acoustic  studies  with  coordinated  sampling.  The  results  of  all  these  efforts 
show  that  ocean  variability  can  considerably  influence  acoustic  propagation  properties  both  in 
time  and  in  space.  Acoustic  propagation  in  shallow  water  is  also  known  to  be  strongly  affected  by 
multiple  seabed  interactions  and  bathymetry  effects  [23].  In  continental  shelf  break  regions,  the 
seabed  and  bathymetry  are  especially  complex  and  not  well  known.  Uncertainties  in  geoacoustic 
modeling  usually  come  from  the  imperfect  knowledge  of  the  true  thicknesses  and  properties 
of  the  sediment  and  rock  layers  in  the  sea  floor.  Studying  and  quantifying  the  sensitivity  of 
our  coupled  predictions  to  key  ocean  and  geoacoustic  uncertainties  is  another  objective  of  our 
present  effort. 

Our  investigations  are  part  of  the  Quantifying,  Predicting,  and  Exploiting  Uncertainty  (QPE) 
initiative  which  aims  to  integrate  probabilistic  performance  prediction,  coupled  ocean-acoustic 
modeling,  multidisciplinary  data  assimilation  and  autonomous  ocean  platforms  to  improve  per¬ 
formance  prediction  and  reduce  detection  uncertainties.  Our  focus  is  on  the  QPE  2008  pilot 
exercise  [24]  carried  out  on  the  continental  shelf  and  slope  of  northeast  Taiwan  in  September, 
2008.  The  main  goals  of  this  2008  exercise  were  to  test  systems  and  methodologies  to  measure 
and  forecast  the  baseline  regional  ocean  variabilities  and  uncertainties,  as  well  as  their  impacts 
on  low-frequency  (100  to  lOOOHz)  acoustic  propagation  conditions.  Two  main  locations  were 
chosen  for  the  exercise,  one  on  the  continental  shelf  and  the  other  on  the  shelfbreak.  Interestingly, 
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the  transmission  loss  over  ranges  of  5  to  6  km  over  eircular  areas  in  these  two  regions  were 
observed  to  be  very  different.  This  oeeurred  even  though  the  two  eireular  areas  overlapped  in 
space.  For  example,  no  signal  was  obtained  when  the  source  was  over  the  deeper  water  on  the 
shelfbreak,  and  the  measured  mean  and  variability  of  the  transmission  loss  differed  greatly  in 
the  two  regions.  In  real-time  and  after  the  experiment,  the  community  could  not  explain  the 
causes  of  these  differences.  To  understand  and  quantify  these  observed  behaviors,  both  acoustic 
and  environmental  dynamical  modeling  are  required. 

In  the  present  study,  we  utilize  coupled  oceanographic(4D)  and  acoustic(Nx2D)  fullfield 
simulations  to  explain  and  quantify  the  mean  and  dynamic  variability  of  mid-frequency  sound 
transmission  losses  observed  during  the  QPE  2008  Pilot  experiment  northeast  of  Taiwan.  To 
do  so,  coupled  model  estimates  are  for  the  first  time  compared  to  the  underwater  acoustic 
observations  available.  We  also  investigate  the  sensitivity  of  the  TL  estimation  to  uncertainties 
in  geo-acoustic  parameters,  bathymetry,  sediment  layer  thicknesses  and  initial  transports  in  the 
Taiwan  Strait  as  well  as  to  background  sound  speed  variability  involving  tidal  effects. 

In  what  follows,  we  first  outline  the  QPE  2008  Pilot  experiment  and  the  corresponding 
observations  and  coupled  modeling  systems  (Sect.  II).  In  Sect.  Ill,  the  transmission  loss  data 
and  its  variability  with  respect  to  positions  and  bearing  angles  are  compared  with  our  coupled 
oceanographic-acoustic  modeling  predictions.  The  uncertainties  in  the  TE  estimates  due  to  var¬ 
ious  factors  are  also  modeled  and  quantified,  and  the  sensitivity  of  these  TE  estimates  to  ocean 
variability  is  computed  and  studied.  The  summary  and  conclusions  are  in  Sect.  IV. 

II.  Real-time  Experiment,  In-situ  Observations  and  Coupled  Modeling  Systems 

The  QPE  2008  Pilot  Cruise  was  carried  out  on  the  continental  shelf  and  slope  of  northeast 
Taiwan  (Eigure  1).  Most  of  the  acoustic  propagation  experiments  were  conducted  inside  the 
rectangular  boxed  region  just  north  of  (25.5°N,122.5°E).  Our  regional  dynamics  and  modeling 
focus  was  the  continental  shelf  and  slope  northeast  of  Taiwan,  and  especially  the  Cold  Dome 
(e.g.  [25]),  its  dynamics,  variabilities  and  uncertainties,  as  well  as  impacts  of  this  environment 
on  low/mid-frequency  (100  to  1000  Hz)  acoustic  propagation.  In  the  region,  a  large  number  of 
ocean  processes  can  occur  simultaneously,  very  energetically  and  on  multiple  scales  [24],  [26]- 
[28].  The  ocean  observations,  the  coupled  ocean-acoustic  modeling  systems  used  for  real-time 
predictions  and  the  acoustic  observations  are  described  next. 
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A.  Ocean  Observations 

The  Pilot  Experiment  involved  three  separate  eruises  to  eolleet  oeean  environmental  data  (see 
Figure  1(a)  for  CTD  stations  positions).  The  solid  dots  are  for  the  ORl  eruises:  Leg  1  was 
operated  between  the  Ilan  Ridge  and  the  Main  Study  Area  during  September  2-4  while  Leg  2 
was  in  the  Main  Study  Area  from  September  6-11.  The  stars  and  eireles  indieate  the  stations 
that  the  vessels  operated  during  August  22-27  (respeetively,  the  OR2  northeast  of  Taiwan  and 
the  OR3  in  the  Taiwan  Strait).  As  soon  as  available,  these  CTD  measurements,  as  well  as  other 
measurements  were  assimilated  in  our  simulations. 

B.  Real-time  Modeling  Systems  and  Coupled  Predictions 

In  real-time,  we  earned  out  sustained  eoupled-oeeanographie(4D)-aeoustie(Nx2D)  fullfield 
simulations  for  a  7-day  long  period  [26],  eaeh  day  seleeting  a  subset  of  these  simulations  as  the 
aeoustie-oeean  foreeast  for  the  next  two  days.  Our  simulations  utilized  the  MIT  Multidiseiplinary 
Simulation,  Estimation  and  Assimilation  System  (MSEAS).  This  system  is  used  for  realistie 
simulations  and  real-time  foreeasts  in  several  regions  of  the  world’s  oeean.  For  the  oeean 
dynamies,  it  employs  a  new  free-surfaee  and  two-way  nested  primitive-equation  eode  [29]  (a 
signifieant  upgrade  of  the  rigid-lid  Harvard  primitive-equation  model  [30],  [31]).  The  MSEAS 
system  also  involves  a  eoastal  objeetive  analysis  seheme  based  on  fast-marehing  methods  [32]; 
an  Optimal  Interpolation  seheme  and  Error  Subspaee  Statistieal  Estimation  system  (ESSE,  [33]- 
[36])  for  data  assimilation,  optimization  and  adaptive  sampling;  new  sehemes  for  uneertainty 
predietions  based  on  dynamieally-orthogonal  equations  [37];  multiple  biologieal  models  [38]; 
and,  several  aeoustie  models. 

The  QPE  oeean  simulations  were  foreed  with  barotropie  tidal  flows  [39]  at  open  oeean 
boundaries  and  with  a  eombination  of  COAMPS  (wind  stress)  and  NOGAPS  (heat-flux,  E-P) 
atmospherie  foreing  at  the  air-sea  interfaee.  Fields  were  initialized  with  the  OR2  and  OR3  CTD 
data  merged  with  a  summer  elimatology  ereated  using  June- August  profiles  and  the  HydroBase2 
software  [40].  The  bathymetry  used  was  the  NCOR  bathymetry  [41].  The  ORl  CTD  and  the 
Seasoar  data  were  assimilated  when  available.  During  the  preparation  of  the  exereise,  we  found 
that  the  net  transport  between  Taiwan  and  mainland  China  ean  have  a  signifieant  influenee  on 
the  formation  and  strength  of  the  eold  dome,  espeeially  just  north  of  Taiwan.  To  aeeount  for  the 
eorresponding  uneertainty  in  real-time,  2-3  different  initial  transport  oases  were  utilized  for  eaeh 
oeean  foreeast  (of  oourse,  boundary  oonditions  vary  in  time  with  a  radiation  eondition  [42],  but 
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the  radiation  condition  is  applied  to  the  departure  from  this  initial  estimate  [43]  which  remains 
fixed).  The  skill  of  the  ocean  predictions  obtained  from  each  of  these  different  initial  boundary 
condition  cases  were  computed  in  real-time  by  daily  comparisons  to  the  available  in  situ  data 
and  SST  [44], 

The  ocean  forecasts  and  their  uncertainties  were  used  as  inputs  to  the  acoustic  simulations.  The 
bathymetry  data  were  provided  by  the  Center  for  Coastal  and  Ocean  Mapping  with  up  to  100  m 
resolution  in  the  main  study  region  [45].  Acoustic  simulations  were  performed  with  the  Coupled 
SACLANTCEN  normal  mode  propagation  loss  model  (C-SNAP)  [46]  using  the  sound-speed 
fields  of  MSEAS  as  parameters.  In  real-time,  transmission  loss  forecasts  were  provided  along 
five  20  km-long  acoustic  propagation  paths  which  had  been  pre-selected  to  explore  the  variability 
in  this  region.  Real-time  results  are  presented  in  [26].  Specifics  on  our  acoustic  modeling  and 
results  are  provided  in  Section  III.  As  a  whole,  the  QPE  scientific  team  successfully  resolved 
the  Cold  Dome,  preformed  data  transfer  in  real  time  and  forecasted  oceanographic  and  acoustic 
propagation  conditions  daily,  including  uncertainties  in  ocean  fields  due  to  uncertainty  in  the 
Taiwan  Strait  transport. 

C.  Acoustic  Propagation  Observations 

The  acoustic  propagation  experiments  were  conducted  during  the  leg  2  of  the  Pilot  Cruise  on 
ORE  The  OMAS  vehicles  were  utilized  as  sound  sources.  The  receivers  included  standard  US 
Navy  sonobouys  and  the  NTU  vertical  line  array  (VEA).  Eour  OMAS  events  were  carried  out, 
two  on  the  shelf  along  relatively  constant  headings  and  depths  and  two  in  circular  tracks  around 
the  shelfbreak  region. 

The  first  two  events  conducted  on  the  shelf  occurred  at  approximately  110  m  to  130  m 
depth.  The  OMAS  moved  at  a  speed  of  5  knots.  Our  numerical  modeling  results  of  these  first 
two  events  were  compared  and  fitted  to  the  transmission  loss  observations  at  different  ranges, 
providing  insights  on  the  values  of  the  geoacoustic  model  parameters.  Eits  were  good  but  results 
are  not  shown  here. 

The  other  two  acoustic  events  were  the  two  circular  OMAS  runs:  Event  A  and  B,  as  shown  in 
Eig.  1(b).  The  northern  and  southern  circles  are  the  nominal  Event  A  and  B  tracks,  respectively. 
The  programmed  radius  of  Event  A  was  set  to  5  km  and  that  of  Event  B  to  6  km.  The  OMAS 
vehicles  and  the  sonobuoys  were  set  to  run  at  61  m  in  depth.  The  northern  circle  was  centered  at 
approximately  1 10  m  in  depth,  and  covered  a  relatively  flat,  slightly  sloping  bottom.  The  southern 
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circle  was  positioned  over  the  steeply  breaking  shelf.  Fig.  ure  2(a)  presents  TL  measurements 
collected  from  the  two  runs:  the  blue  curve  denotes  the  1100  Hz  TL  data  for  Event  A,  the  red 
curve  is  the  900  Hz  TL  data  for  Event  B.  Note  the  TL  data  here  is  the  peak  of  the  matched  filer 
output  and  small  linear  range  corrections  are  applied  to  obtain  the  TL  estimation  at  the  range 
of  6  km.  The  TL  of  Event  A  varies  somewhat  with  the  bearing  angle,  showing  an  increased 
TL  in  the  north-eastern  direction.  Interestingly,  during  Event  B  (southern  circle  positioned  over 
the  shelfbreak),  no  signal  was  observed  when  the  OMAS  source  was  in  the  deeper  waters.  In 
addition,  considering  the  arcs  over  the  shelf  only,  the  mean  of  the  TL  data  was  2-5  dB  larger  for 
Event  B  (shelfbreak  run)  than  for  the  arc  of  Event  A  (shelf  run).  The  goals  of  our  modeling  are 
to  explain  and  quantify  such  observations,  and  to  investigate  the  variation  of  the  TL  estimates 
due  to  uncertainties  in  the  sediments,  bathymetry,  and  ocean  water  column  variation. 

III.  Results 

Our  simulations  start  with  the  range  dependent  ocean  sound  speed  fields  provided  by  the 
forecasts  of  our  MIT  MSEAS  system  for  the  period  of  September  I-I2,  2008  (Section  II-B). 
Acoustic  TL  fields  were  computed  using  C-SNAP  in  180  bearing  directions  (2-degree  resolution) 
for  both  Events  A  and  B,  at  frequencies  of  1100  Hz  and  900  Hz,  respectively.  In  these  Nx2D 
computations,  the  two  simulated  OMAS  circles  were  centered  at  (25°N  39.8’  122°E  35.76’) 
and  (25°N  44.55’  122°E  33.90’),  with  radius  of  6  km,  respectively.  These  locations  are  mean 
estimates  of  the  positions  of  the  slightly  moving  receivers’  tracks  [24].  The  sources  and  receivers 
were  set  at  61  m  depth  below  the  sea  surface  to  fit  the  sea  configuration.  Transmission  losses 
were  computed  at  positions  that  were  ±  600  m  in  range  and  ±5  m  around  the  depth  of  the 
receiver.  The  average  predicted  TL  within  these  horizontal  and  depth  ranges  was  then  set  to 
be  our  TL  estimate  because  it  is  a  good  approximation  [47]  to  the  10  percent  bandwidth  of 
frequency  averaging  that  occurred  in  the  real  in  situ  broadband  propagation.  All  of  the  Nx2D 
TL  estimates  that  we  show  in  this  manuscript  are  broadband  estimates,  but  they  are  often  simply 
referred  to  as  TL.  The  first  results  of  our  coupled  ocean-acoustic  simulations  are  illustrated  in 
Ligure  2(b),  with  Event  A  at  12:00,  September  8,  2008,  and  Event  B  at  15:00,  September  8, 
2008  (all  UTC  times).  Shown  are  our  broadband  TL  estimates  for  the  time  and  location  of  the 
two  events,  for  all  bearing  angles.  If  we  consider  only  the  range  of  bearing  angles  for  which 
a  signal  was  received  at  sea  from  the  OMAS  runs,  our  mean  TL  estimates  for  Event  A  and 
Event  B  are  about  77-81  dB  and  79-85  dB,  respectively.  These  predicted  numbers  are  in  good 
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agreement  with  the  mean  data  received,  as  illustrated  on  Fig.  2(a).  In  fact,  this  is  the  first  time 
that  coupled  ocean-acoustic  predictions  are  validated  using  in  situ  acoustic  data.  The  measured 
and  predicted  mean  TL  are  within  1  to  3dB,  which  is  within  the  uncertainties  of  our  coupled 
acoustic-ocean  simulations  (Sect.  III).  In  addition  to  the  mean  values,  spatial  patterns  are  also 
captured.  For  example,  for  Event  A,  one  notes  the  increased  TL  in  the  north-eastern  direction 
(225°  to  45°),  in  accord  with  the  data  (Fig.  2(a)):  this  is  linked  to  ocean  variability  as  shown  in 
Sect.  III-D.  For  Event  B,  our  simulated  mean  TL  on  the  shelf  is  also  lower  by  2-5  dB  than  that 
of  Event  A  along  the  same  angles,  as  was  observed:  this  result  is  linked  to  seabed  properties 
(Sect.  III-A). 

The  bathymetry  data  (see  Eig.l(c))  along  bearing  angles  of  120  to  230°  shows  depth  variations 
from  150  m  to  more  than  300  m.  In  our  coupled  modeling  results,  this  leads  to  a  mean  TL 
estimate  for  Event  B  that  is  higher  by  15-20  dB  than  the  mean  TL  of  Event  A  along  the  same 
angles  (Eigure  2(b)).  We  expect  that  this  explains  why  there  was  no  reception  during  the  Event 
B  run  when  the  OMAS  went  over  this  deeper  shelfbreak  region.  At  these  angles,  this  leads  to 
an  up-slope  transmission  which  causes  extensive  transmission  loss.  However,  since  the  ocean, 
seabed  and  bathymetry  are  all  uncertain,  it  is  also  possible  that  some  water  column,  seabed  or 
bathymetric  features  were  the  causes,  or  were  at  least  in  part  responsible,  of  this  behavior.  To 
study  and  quantify  these  possibilities,  we  investigate  next  the  sensitivity  of  our  broadband  TL 
estimates  to  uncertainties  in  geo-acoustic  parameters  (sediment  thickness  and  sediment  types), 
bathymetry,  and  variability  of  the  ocean  sound  speed  fields  (Taiwan  strait  transport  and  tidal 
dynamics).  We  also  quantify  the  effects  of  these  uncertainties  on  the  mean  TL  estimates.  Eor 
each  uncertainty  factor,  we  compute  the  standard  deviations  of  our  TL  estimates  (assuming  all 
other  factors  constant  and  perfectly  known).  Of  course,  these  standard  deviations  differ  from  the 
standard  deviations  of  the  TL  measured  at  sea.  This  is  because  the  standard  deviations  of  the 
measured  TL  correspond  only  to  variations  in  the  ocean  and  in  the  positions/bearing  during  the 
given  OMAS  circle  run  (assuming  the  bathy  and  seabed  properties  are  fixed  during  that  time). 

A.  TL  sensitivities  to  uncertainties  in  the  Geoacoustic  model  -  sediment  types  and  thicknesses 

Uncertain  factors  in  the  geoacoustic  model  that  lead  to  uncertainties  in  the  TL  estimates  include 
the  thicknesses  and  properties  of  the  sediment  and  rock  layers  in  the  sea  floor.  At  higher  sound 
frequencies,  the  first  few  to  tens  of  meters  of  sediments  matter  most.  At  lower  frequencies,  the 
whole  sediment  column  and  underlying  rocks  often  matter.  A  challenge  for  precise  modeling  of 
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sound  propagation  in  this  region  is  the  laek  of  knowledge  of  the  geo-aeoustie  parameters.  For  our 
frequencies  (900  and  1110  Hz),  we  evaluated  and  compared  more  than  thirty  sediment  models. 
Here  we  only  report  a  few  of  them  including  the  properties  of  sand,  muddy-sand,  silty-clay  and 
hybrid  models  in  Tables  I  and  II.  Specifically,  Table  I  gives  three  sediment  types  including  an 
additional  half-space  basement,  while  Table  II  gives  our  best-fit  so  far,  an  hybrid  and  depth- 
dependent  model.  For  each  sediment,  note  that  the  first  line  corresponds  to  the  optimum  values 
we  obtained  using  the  hybrid  model,  while  the  second  line  indicates  the  deviation  ranges  of  the 
parameter  values  that  we  tested  on.  The  first  three  sediments  followed  Hamilton’s  models  [48], 
[49].  The  hybrid  depth-dependent  models  were  the  best- fit  results  of  our  multiple  comparisons 
to  data  aiming  to  optimize  parameters  (three  sediment  types  and  two  transition  depths)  and  of 
interactions  with  Heaney,  Holland  and  Chen  (personal  communications).  A  similar  hybrid  model 
was  also  obtained  by  Heaney  [50]. 

In  Figure  3,  we  illustrate  how  sediment  layers  affect  our  TL  estimates  for  both  Events  A  and 
B,  showing  the  results  for  four  of  the  models  we  evaluated  (see  Tables  I  and  II).  The  sediment 
thickness  is  set  to  20  m  in  all  cases.  We  find  that  the  four  models  have  different  impacts  on  the 
shelf  and  shelfbreak. 

For  Event  A  on  the  shelf,  the  four  sediments  lead  to  distributions  of  TE  variation  with  bearing 
angles  that  are  very  similar.  The  sediment  types  shown  lead  to  overall  isotropic  TE  differences 
of  5  to  15  dB  in  all  directions.  The  average  standard  deviation  is  9.55  dB. 

Eor  Event  B  on  the  shelfbreak,  the  average  standard  deviation  is  smaller,  5.81  dB,  than  for 
Event  A.  The  shapes  of  the  four  TE  estimates  are  still  similar,  but  a  bit  less  than  for  Event  A.  The 
TE  deviations  due  to  different  sediments  are  estimated  to  be  lesser  along  the  shelfbreak  (e.g.  5  dB 
or  less)  but  a  bit  larger  towards  the  shelf  (e.g.  about  7-8  dB  along  330°)  and  for  bearing  angles 
of  120°  to  240°  which  is  over  the  deeper  bathymetry  region.  In  fact,  the  steepest  bathymetries 
lead  to  a  very  localized  amplification  of  TE  uncertainties  due  to  sediment  uncertainties:  TE 
deviations  are  maximum  (30  dB)  along  the  bearing  angle  of  180°.  Nonetheless,  all  TE  fields  of 
Event  B  for  bearing  angles  of  120°  to  240°  clearly  show  a  much  higher  TE  than  in  Event  A  (on 
average,  by  about  20  db),  as  was  measured  at  sea  (see  Eig.  2(a)).  This  implies  that  the  sediment 
properties  cannot  alone  explain  the  observed  higher  loss. 

We  also  investigated  the  effects  that  uncertainties  in  the  thicknesses  of  sediment  layers  have 
on  our  TE  estimates.  The  results  (plots  not  shown  here)  indicate  that  variations  of  5  to  25  m 
in  the  sediment  thicknesses  have  some  limited  effects  only  on  the  up-slope  transmissions  for 
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the  circle  runs  of  both  Event  A  (less  than  2  dB  uncertainty)  and  Event  B  (less  than  2-3  dB 
uncertainty,  excepted  in  the  deeper  region  where  it  reaches  3-6  dB). 

An  important  result  of  these  sensitivity  studies  to  uncertainties  in  both  sediment  types  and 
thicknesses  is  that  neither  of  them  can  alone  explain  the  lack  of  transmission  for  Event  B  when 
the  OMAS  source  moved  over  the  deeper  shelfbreak.  A  combination  of  both  uncertainties  (not 
shown)  still  does  not  change  this  conclusion.  Therefore,  the  bathymetry  itself  is  found  responsible 
for  the  lack  of  transmission.  A  second  result  is  that  although  the  full  geo-acoustic  inversion  has 
not  yet  been  completed  for  this  site,  our  simple  optimized  hybrid  sediment  model  and  sediment 
thickness  (20  m)  with  our  4D  simulations  of  the  ocean  sound  speed  lead  to  predicted  estimates 
of  the  broadband  TE  that  are  in  good  agreement  with  the  measurements  (see  Eig.  2(b)). 

B.  TL  sensitivities  to  uncertainties  in  bathymetry 

Since  we  found  that  bathymetry  is  likely  responsible  for  causing  the  loss  of  acoustic  reception 
in  Event  B  over  the  shelfbreak,  we  estimated  the  TE  variations  due  to  uncertainties  in  the 
bathymetry.  By  comparing  bathymetries  of  different  resolutions  and  different  origins,  we  obtained 
a  general  characteristic:  errors  in  bathymetry  are  a  function  of  the  depth  and  slope.  Therefore, 
we  derived  a  simple  model  of  the  probability  density  of  the  baythmetry  around  its  best  estimate 
(100m  resolution).  This  model  consists  of  a  zero  mean  Gaussian  noise  with  variance  a  function 
of  the  slope  and  depth,  i.e., 

D{r-,uj)  =  D{r)[l  +  a{r)  ■  u],  (1) 

where  D{r;u)  is  a  stochastic  process,  12  is  the  sample  space  containing  the  set  of  elementary 
events  w  G  12  ~  A^(0, 1),  r  is  range  and  a{r)  =  e-S{r)  where  the  normalized  slope  S{r)  C  [0  1]. 
The  simulation  results  with  different  realizations  of  this  bathymetric  uncertainty  model  are  shown 
in  Eigure  4.  Each  TE  circular  curve  corresponds  to  a  bathymetric  noise  with  different  percentages 
of  local  depth  and  slope. 

As  a  whole,  for  these  levels  of  bathymetric  uncertainties,  we  find  that  the  uncertainties  of 
our  TE  estimates  are  proportional  to  uncertainties  in  the  bathymetry.  Another  result  is  that 
the  general  shape  of  the  TE  circular  curves  with  bearing  angles  remains  as  obtained  before. 
Therefore,  bathymetric  uncertainties  cannot  explain  the  observed  increase  in  TE  over  the  deeper 
ocean  region.  Again,  bathymetry  itself  provides  the  explanation. 
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C.  TL  sensitivities  to  uncertainties  in  the  ocean  water  column 

In  the  following,  we  analyze  effeets  of  uneertainties  in  the  mean  water  eolumn  on  our  aeoustie 
transmission  loss  predietions.  We  foeus  on  uneertainties  in  the  transport  through  the  Taiwan  Strait 
(oeean  between  mainland  China  and  Taiwan,  e.g.  [51]).  Even  though  the  operational  aeoustie 
area  is  more  than  100  km  away  northeast  of  the  Taiwan  Straits  eenterline,  we  had  diseovered 
(see  Seet.  II-B)  that  the  strength  of  the  Straits  transport  eould  have  a  signifieant  influenee  on 
the  Cold  Dome  region,  espeeially  just  north  of  Taiwan  [26].  If  the  initial  Taiwan  Strait  transport 
was  northward  and  of  the  order  of  ISv  or  more,  our  simulations  showed  that  the  Kuroshio  had 
limited  intrusion  over  the  shelfbreak  northeast  of  Taiwan  and  the  Cold  Dome  did  not  really 
form.  However,  if  the  initial  Taiwan  Strait  transport  was  weaker  or  southward,  we  found  that 
deeper  and  eooler  Kuroshio  waters  upwelled  on  the  shelf  and  the  Cold  Dome  eould  form,  the 
strength  of  this  formation  being  a  funetion  of  the  transport.  For  the  aeou  sties  in  the  operational 
area  northeast  of  Taiwan  (see  Fig.l),  these  different  oeean  foreing  ehange  the  mean  sound  speed 
field.  They  ean  thus  have  signifieant  impaets  on  the  TF  estimates. 

In  real-time,  we  eonsidered  three  different  initial  transport  eonditions  through  the  Taiwan 
Strait:  an  initial  transport  of  1  Sv  northward,  of  1  Sv  southward,  and  of  0  Sv.  Figure  5  shows 
sound  speed  fields  in  the  QPE-Pilot  experiment  region  at  30  m,  60  m  and  90  m  depth  (lines  1, 
2  and  3  on  Fig.  5,  respeetively)  with  different  initial  transport  eonditions  through  the  Taiwan 
Strait  (eolumns  1,  2  and  3  on  Fig.  5,  respeetively).  The  main  aeoustie  researeh  area  is  denoted 
by  a  white  reetangle  frame.  The  differenees  of  sound  speed  at  the  depth  of  30  m,  and  60  m 
depth  are  easily  notieeable,  but  not  as  apparent  at  the  depth  of  90  m  in  part  beeause  the  range 
of  sound-speeds  present  in  the  region  is  larger  at  this  depth. 

Figure  6  illustrates  our  Nx2D  aeoustie  simulations  for  Events  A  and  B,  using  the  three  different 
sound  speed  baekgrounds  illustrated  by  the  eolumns  of  Fig.  5.  The  geo-aeoustie  model  in  all 
eases  was  the  hybrid  depth-dependent  model  with  a  20-m  thickness  of  sediment  layer  (see 
Sect.  Ill- A).  The  three  4D  sound  speed  fields  were  predicted  by  our  MSEAS  ocean  forecast 
system:  each  correspond  to  a  different  initial  transport  condition  between  Taiwan  and  mainland 
China  (±1  Sv  and  0  Sv).  The  prediction  skill  for  each  of  these  initial  conditions  was  computed 
in  real-time  daily,  based  on  comparisons  with  the  available  in  situ  data  and  SST  [44]  for  that  day. 
On  September  8,  2008,  the  1  Sv  South  initial  transport  condition  was  chosen  as  the  best  ocean 
forecast.  We  find  that  the  three  ocean  fields  cause  different  effects  on  the  acoustic  transmissions 
of  Events  A  and  B.  On  Sep  8th  (the  day  when  Events  were  run  at  sea),  the  uncertainty  in  the 
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ocean  state  due  to  the  uncertainty  in  the  initial  eonditions  of  the  Taiwan  Strait  transport  has  more 
impacts  on  our  broadband  TL  estimates  for  the  shelf  break  Event  B  than  for  the  shelf  Event  A. 
Eor  the  shelf  Event  A,  our  estimated  TE  variations  along  the  different  bearing  angles  reached 
about  1  to  3  dB,  with  a  bearing-averaged  standard  deviation  of  0.58dB.  However,  for  the  shelf 
break  Event  B,  our  TE  estimates  varied  almost  evenly  in  all  direetions  up  to  about  2  to  5  dB, 
with  a  bearing-averaged  standard  deviation  of  1.21dB. 

In  eonelusion,  these  results  first  indieate  that  uneertainties  in  the  Taiwan  Strait  transport  lead  to 
broadband  TE  uneertainties  that  are  smaller  on  the  shelf  than  on  the  shelfbreak.  This  is  beeause 
the  shelfbreak  ocean  dynamies  is  more  sensitive  to  this  transport  than  the  shelf  dynamies.  Another 
result  is  that  the  shape  of  our  broadband  TE  eurve  does  not  vary  much  with  these  different  oeean 
foreeasts:  it  is  mainly  the  mean  TE  whieh  is  affeeted.  This  indieates  that  the  laek  of  signal 
reeeived  for  Event  B  when  the  OMAS  passed  over  the  deeper  shelfbreak  region  is  not  due  to 
these  oeean  uncertainties  (themselves  due  to  the  uneertain  initial  eonditions  of  the  Taiwan  Strait 
transport).  The  bathymetry  remains  the  explanation. 

D.  TL  sensitivities  to  tidal  dynamics  in  the  ocean  water  column 

We  now  examine  the  variations  of  the  transmission  loss  field  due  to  variability  in  the  oeean, 
speeifically  tidal  effeets.  We  foeus  on  internal  tidal  variability  beeause  our  data-assimilative  sim¬ 
ulations  revealed  that  semi-diurnal  frequeneies  in  the  oeean  interior  led  to  signifieant  variations 
in  the  sound  energy  transmitted.  Eirst,  we  charaeterize  this  oeean  variability,  then  we  study  its 
impaets  on  the  broadband  sound  transmission. 

Our  simulations  of  the  sound-speed  field  in  the  region  over  a  24-hour  period  are  illustrated 
on  Eigure  7,  for  the  day  on  whieh  the  two  acoustie  Events  were  measured  at  sea.  Shown  is 
the  sound-speed  variability  at  30  m  depth  (in  the  upper  layers  of  the  main  thermoeline),  every 
3  hours  from  00:00Z  Sep.  8  to  00:00Z  Sep.  9,  2008.  Panel  (a)  is  the  total  sound-speed  map 
on  00:00Z  Sep.  8  while  the  other  8  Panels  (b)-(i)  are  the  differenees  between  the  subsequent 
3-hourly  maps  and  this  00:00Z  Sep.  8  map.  Panels  are  ordered  aeeording  to  time  sueh  that:  those 
on  the  main  diagonal  (a,  e  and  i)  are  12  hours  apart  from  each  other  (elose  to  be  in  phase  for 
the  semi-diurnal  tidal  period);  those  on  the  first  off-diagonals  (b  and  f,  d  and  h)  are  3-hours  off 
from  those  on  the  main  diagonal;  and,  those  on  the  seeond  off-diagonals  (e  and  g)  are  6-hours 
off  from  those  on  the  main  diagonal  (elose  to  be  in  opposition  of  phase  for  the  semi-diurnal 
period). 
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The  predicted  sound-speed  variability  maps  clearly  show  the  presence  of  internal  tides  ori¬ 
ented  along  the  shelfbreak  over  more  than  a  200  km  extent.  They  propagate  up-slope  and  our 
simulations  predict  that  their  wavelength  is  around  40  to  80  km.  These  maps  also  indicate  that 
the  shelfbreak  internal  tide  patterns  affecting  our  acoustic  region  (the  white  rectangle)  are  weaker 
and  not  aligned  with  the  tide  patterns  around  the  tip  of  Taiwan.  Other  internal  patterns  are  also 
present  in  the  Taiwan  strait.  At  other  depths  (e.g.  60  m,  90  m,  not  shown),  variations  in  time 
reveal  similar  properties.  Overall,  we  find  that  the  amplitude  of  shelfbreak  internal  tides  decays 
as  the  waves  propagate  up  the  shelf  but  remain  significant  up  to  about  the  80  m  isobath.  After 
24  hours  (Panel  i),  predicted  variations  also  show  effects  of  mesoscale  features  (eddies,  fronts, 
etc.). 

Figures  8  and  9  illustrate  the  sound  speed  variability  in  vertical  sections  along  and  across 
the  shelfbreak  direction,  respectively.  The  locations  of  the  two  sections  and  the  acoustic  circles 
of  Events  A  and  B  are  plotted  inside  the  white  rectangular  box  drawn  on  the  maps  of  Figs.  5 
and  7.  Panels  are  still  shown  every  3  hours  during  Sep.  8  (UTC).  Panel  (a)  is  the  total  sound- 
speed  section  on  00:00Z  Sep.  8  while  the  other  8  Panels  (b)-(i)  are  the  differences  between  the 
subsequent  3-hourly  sound-speeds  and  this  00:00Z  Sep.  8  sound-speed. 

Within  the  section  along  the  shelfbreak  (Fig.  8),  one  detects  the  crest  and  trough  of  wave 
patterns:  they  are  aligned  with  the  section  and  concentrated  around  the  main  thermocline.  As 
shown  by  the  3-hourly  difference  fields  (Panels  b-i),  the  patterns  are  clearly  depth-dependent 
(baroclinic)  and  of  tidal  period  (e.g.  see  the  negative  and  positive  anomalies  on  the  upper  and 
lower  first-diagonals  of  Fig.  8,  respectively).  This  confirms  the  internal  tide  pattern.  This  signal 
is  not  small,  the  largest  amplitudes  are  around  4  m/s.  The  vertical  extent  is  from  the  bottom 
of  the  mixed  layer  (20  m  below  the  ocean  free  surface)  to  about  10  m  above  the  seafloor.  The 
isotherms  in  the  thermocline  oscillate  up  and  down  with  an  amplitude  of  about  lO-to-20  m 
(note  that  vertical  oscillations  of  different  isotherms  are  not  exactly  in  phase  as  indicated  by  the 
baroclinic  behavior). 

Within  the  section  across  the  shelfbreak  (Fig.  9),  our  simulations  show  the  train  of  wave 
patterns  moving  up  the  slope  and  shelf,  across  the  shelfbreak.  Again,  the  diagonal  Panels  are 
close  to  be  in  phase,  while  the  most  (second)  off-diagonal  Panels  are  close  to  be  in  opposition 
of  phase  with  these  diagonal  Panels  for  the  semi-diurnal  period.  The  internal  tide  wavelength  at 
this  location  and  time  is  found  to  be  within  40  to  60  km.  The  largest  amplitudes  of  the  waves 
in  this  section  are  around  4-5  m/s. 
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The  internal  tide  variability  that  our  simulations  revealed  (Figs.  8-9)  was  found  to  have 
signifieant  impaets  on  our  simulated  broadband  sound  transmission  loss  in  the  aeoustie  region. 
These  results  are  illustrated  now.  Speeifieally,  the  variations  in  time  of  our  Nx2D  TL  simulations 
of  Events  A  and  B  are  shown  in  Fig.  10  again  for  the  day  (Sep.  8th,  2008)  when  the  two  OMAS 
runs  oeeurred  at  sea.  The  nine  different  times  still  refer  to  3-hour  intervals.  Panel  (a)  is  the 
total  broadband  TL  on  00:00  Sep.  8  (from  65  to  105  dB,  w.r.t.  aim  souree).  The  other  8 
Panels  (b)-(i)  are  the  differenees  between  the  subsequent  3-hourly  TLs  and  this  00:00  Sep.  8 
TL.  Note  that  the  scale  for  these  latter  differences  is  -10  dB  to  10  dB. 

Concentrating  first  on  the  shelf  (Event  A,  blue  curves),  we  predict  relatively  isotropic  broad¬ 
band  TL  variations  (much  more  isotropic  than  for  Event  B,  red  curves):  in  other  words,  thermo- 
cline  motions  due  to  internal  tides  lead  to  TL  variations  on  the  shelf  that  are  more  uniform  with 
bearing  angles  than  they  are  across  the  shelfbreak  (Event  B).  This  is  because  the  6  km  range  of 
OMAS  transmissions  are  smaller  than  half-a-wavelength  of  the  internal  tides  (see  Figs.  8  and  9). 
Hence,  for  the  relatively  flat  shelf  bathymetry,  the  oscillating  depth  of  the  thermocline  governs  the 
TL  over  6  km.  On  the  shelf,  the  anisotropy  of  our  predicted  TL  over  6  km  range  arises  because 
of  the  recurring  slopes  of  the  thermocline:  see  the  similarity  between  the  differences  of  the  first 
lower-diagonal  (Panels  d  and  h),  first  upper-diagonal  (Panels  b  and  f)  and  second  off-diagonal 
(Panels  c  and  g).  Note  also  that  differences  are  close  to  zero  on  the  main  diagonal  (Panels  e 
and  i:  on  bearing-average,  -0.5  and  0.9dB,  respectively).  On  bearing-average,  the  differences  are 
largest  for  Panels  d  and  h:  3.1  and  4.6dB,  respectively.  TL  variability  and  sound-speed  variability 
are  not  exactly  in  phase,  in  part  due  to  the  nonlinear  coupling  and  to  the  4D  processes. 

Over  the  shelfbreak  (Event  B,  red  curves),  the  differences  between  TL  curves  at  different 
times  are  anisotropic  and  rapidly  variable  with  bearing  angles.  At  sea,  it  would  thus  be  harder 
to  distinguish  coherent  internal  tide  effects  above  the  shelbreak.  This  is  because  the  steep 
and  complex  bathymetry  amplifies  the  internal  tide  oscillations  into  intricate  bearing-dependent 
patterns:  bottom  bounces  at  slightly  different  locations  actually  occur  at  very  different  depths 
and  slope  angles,  leading  to  very  different  transmissions.  This  leads  to  standard  deviations  that 
are  larger  over  the  shelfbreak  (Event  B)  than  over  the  shelf  (Event  A),  e.g.  3.29  dB  vs.  2.21  dB. 
However,  the  bearing-averaged  differences  are  much  smaller  over  the  shelfbreak  than  over  the 
shelf,  as  shown  above.  Over  the  shelbreak,  the  largest  bearing-averaged  differences  are  only 
-1.4  dB  and  1.4  dB,  for  Panels  (c)  and  (i),  respectively,  while  on  the  shelf,  they  reach  4.6  dB  for 
Panel  (h).  Our  simulation  results  are  overall  in  accord  with  observations  of  [22]  for  the  South 
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China  Sea  and  canonical  simulations  of  [52],  even  though  amplitudes  are  different. 

Finally,  one  of  our  main  results  is  not  affected  by  these  TL  variations  (Panels  b-i)  due  to 
internal  tides:  the  mean  TL  over  the  shelfbreak  (Event  B)  is  much  higher  between  120  and  260 
degrees  than  the  mean  TL  over  the  shelf  (Event  A),  regardless  of  the  internal  tide  properties 
(phase,  etc).  This  confirms  that  the  steep  bathymetry  (and  not  the  internal  tides)  explains  the 
lack  of  transmissions  within  these  angles  at  the  location  and  time  of  this  shelfbreak  Event  B. 

To  show  detailed  effects  of  internal  tides  on  the  TL,  full  vertical  sections  are  plotted  on 
Eigure  11.  They  correspond  to  the  0"  due  north  section  in  Eig.  10.  Specifically,  they  are:  a)  and 
b)  sound  speed  sections  at  12:00  and  21:00  on  Sep.  8th;  c)  difference  between  these  two  sections 
(9  h  apart);  d)  and  e)  full  field  TL  sections  for  these  two  times;  and  f),  same  TL  estimates,  but 
at  a  receiver  depth  of  61  meter.  As  seen  from  Panels  (a-c),  after  nine  hours,  the  thermocline  is 
up  by  more  than  10  m  (note  Panel  c  is  close  to  be  a  zoom  in  Pig.  9d,  around  120  m  depth). 
The  vertical  motion  of  the  thermocline  modifies  the  sound-speed  interface  structure.  As  seen  on 
Panels  (d-e),  this  increases  the  width  of  the  TL  convergence  zone  from  1.6  km  to  2  km  and 
causes  more  than  10  dB  variation  in  TL  at  61  m  and  at  6  km  range,  as  shown  in  Pig.  11(f). 

On  Pigure  12,  we  show,  for  both  Events  A  and  B,  our  predicted  bearing-averaged  TL  as  a 
function  of  time  during  Sep.  8th  (UTC).  Interestingly,  we  find  that  the  TL  around  the  times  of 
Event  A  (12  UTC)  and,  to  a  lesser  extent.  Event  B  (15  UTC),  are  within  an  internal  tide  period 
corresponding  to  higher  loss,  see  also  Pig.  10.  Lor  Event  A,  TL  variations  are  due  mostly  to 
thermocline  oscillations  driven  by  internal  tides  and  amplitudes  are  large,  reaching  5  dB.  This 
is  significant  because  our  TL  is  a  broadband  estimate  computed  by  range  and  depth  averaging 
(see  start  of  Sect.  III).  Lor  the  single  frequency  TL  (continuous  wave  -  cw,  computed  prior  to 
averaging),  we  find  that  these  variations  have  amplitudes  of  20  dB.  The  averaging  reduces  the 
amplitudes,  but  they  are  still  5  dB.  In  conclusion,  effects  of  internal  tides  on  the  shelf  correspond 
to  the  largest  TL  sensitivity  that  we  have  found  among  all  factors  studied  above,  almost  in  par 
with  those  due  to  very  different  seabed  properties  (see  Pig.  3). 

IV.  Summary  and  Conclusions 

This  study  was  motivated  by  the  striking  difference  in  acoustic  transmission  data  collected  on 
the  shelf  and  shelf  break  in  the  north-eastern  Taiwan  region  within  the  context  of  the  Quantifying, 
Predicting  and  Exploiting  (QPE)  uncertainty  Pilot  Experiment  [24].  On  the  shelf,  the  mean 
acoustic  transmission  from  a  sound  source  moving  for  6  hours  along  a  circular  track  showed 
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little  TL  variation  with  respect  to  bearing  angles.  The  largest  variations  there  were  increased  loss 
along  the  northeastern  direction  but  with  small  amplitudes.  However,  on  the  shelfbreak,  the  data 
for  the  same  type  of  circular  tracks  showed  no  transmission  when  the  moving  source  was  in  the 
deeper  waters  [53].  In  this  study,  we  explained,  quantified,  and  modeled  the  dynamical  causes 
of  these  observations.  To  do  so,  coupled  oceanographic(4D)-acoustic(Nx2D)  field  estimation 
with  ocean  data  assimilation  was  employed  and,  for  the  first  time,  compared  to  the  observed 
acoustic  data  available.  Predictive  skill,  uncertainties  and  variability  in  time  and  space  were 
all  quantified.  Specifically,  we  studied  the  sensitivity  of  our  results  to  uncertainties  in  a  varied 
set  of  factors,  including  geo-acoustic  parameters,  bottom  layer  thickness,  bathymetry  and  initial 
transport  conditions  in  the  region.  We  also  quantified  spatial  and  temporal  effects  of  internal  tide 
forcing  on  the  the  sound  speed  field  and  acoustic  transmission  field,  both  on  the  shelf  and  over 
the  shelfbreak. 

From  our  coupled  simulations,  we  revealed  that  there  were  about  20  dB  higher  transmission 
losses  in  Event  B  when  the  sound  source  was  in  the  deep  water  region  on  the  shelfbreak. 
We  showed  that  this  dramatic  increase  in  transmission  loss  is  due  to  the  deeper  and  complex 
bathymetry:  our  uncertainty  and  sensitivity  studies  ruled  out  all  other  factors.  This  is  because 
when  the  sound  source  moved  to  the  steep  shelfbreak  region  with  up-slope  transmission,  the 
sound  propagated  down  at  depth  trapping  much  energy  into  the  bottom  and  resulting  in  a  TL 
much  higher  than  the  TL  in  the  reverse  direction  with  down-slope  propagation  over  the  shallower 
shelf. 

Secondly,  we  investigated  uncertainties  in  geo-acoustic  parameters  in  both  the  shelf  and 
shelfbreak  acoustic  simulations.  We  found  that  the  sediment  layer’s  properties  led  to  larger 
but  isotropic  variations  on  the  shelf  (Event  A)  and  smaller  but  more  anisotropic  variations  over 
the  shelfbreak  (Event  B).  We  also  showed  that  our  hybrid  depth-dependent  geo-acoustic  model 
for  the  shelfbreak  and  sand  model  for  the  shelf  better  agreed  with  the  observations  than  other 
models.  Our  investigations  on  the  thickness  of  sediments  showed  that  it  had  some  limited  effects 
(up  to  about  2  dB),  but  only  on  the  up-slope  transmissions  on  both  Events  A  and  B. 

By  comparing  different  sources  and  resolutions  of  accurate  bathymetric  data  sets,  we  derived 
a  statistical  model  of  bathymetric  uncertainties.  Our  resulting  sensitivity  studies  on  these  un¬ 
certainties  revealed  that  TL  uncertainties  (up  to  1  dB  std)  were  proportional  to  the  1  to  4% 
bathymetric  uncertainties  but  that  the  mean  TL  curves  did  not  change  their  shape  as  a  function 
of  bearing  angles. 
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Our  coupled  oceanographic-acoustic  modeling  studies  uncovered  a  surprising  result:  initial 
transport  conditions  in  the  Taiwan  Strait  can  affect  acoustic  transmissions  downstream  more 
than  100  km  away.  We  found  that  it  affected  the  sheflbreak  region  (1  db  std,  up  to  5db)  more 
than  the  shelf  region  (0.5db  std,  up  to  3dB).  This  is  because  the  shelfbreak  ocean  dynamics  is 
more  sensitive  to  this  transport  than  the  shelf  dynamics. 

Using  our  data-assimilative  ocean- acoustic  full  field  modeling  approach,  we  also  described  and 
studied  internal  tide  patterns  over  the  shelf  and  slope  region  north  of  Taiwan  and  quantified  their 
effects  on  transmission  loss  in  the  acoustic  region.  The  wave  patterns  revealed  were  clearly  depth- 
dependent  (baroclinic)  and  close  to  tidal  period,  affecting  the  whole  thermocline  and  extending 
over  more  than  200  km  on  the  shelf.  One  type  of  waves  found  were  shelfbreak  internal  tides 
propagating  up-slope  with  wavelengths  around  40  to  80  km.  Isotherms  were  simulated  to  oscillate 
up  and  down  with  amplitudes  of  lO-to-20  m.  We  discovered  that  these  variations  had  significant 
effects  on  the  transmission  losses  over  5-6  km  range:  on  the  shelf,  variations  of  broadband  TL 
estimates  were  more  isotropic  and  relatively  large  (up  to  5  dB)  than  on  the  shelfbreak  where 
variations  varied  rapidly  with  bearing  angles  due  to  the  steeper  and  more  complex  bathymetry. 
Overall,  internal  tides  on  the  shelf  led  to  the  largest  TL  sensitivity  that  we  simulated.  For  ranges 
of  0(1-10  km),  this  is  a  signal  that  can  be  exploited. 

There  are  several  possible  future  research  activities  linked  to  our  results.  For  example,  we 
have  not  discussed  the  effects  of  the  drifts  of  the  OMAS  sound  source  and  of  the  sonobuoy 
receiver  during  the  experiment  since  these  drifts  were  mostly  accounted  for  in  the  processed 
data  and  in  our  simulations  set-up.  For  the  up-slope  transmissions,  we  also  did  not  discuss  back 
scattering  effects  that  could  arise  due  to  possibly  large-rock  formations  and  volcanic  effects  over 
the  shelfbreak.  The  azimuth  angle  coupling  of  3D  sound  propagation  would  also  affect  our  TL 
estimates,  especially  further  near  the  Canyon  regions.  Our  coupled  results  helped  prepare  for  the 
QPE  intensive  observation  period  in  the  region.  In  general,  we  expect  that  they  will  be  useful 
for  diverse  applications,  including  determining  optimum  sampling  plans,  exploiting  features  and 
improving  4D  coupled  acoustic-ocean  data  assimilation  and  predictions. 
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Cp  (m/s) 

Rho  (g/cm®) 

Attn  (dB/A) 

Sand 

1562 

1.9 

0.90 

1550~1650 

0.25~  1.0 

Muddy-Sand 

1549 

1.488 

1.15 

1540~1575 

0.4  ~1.25 

Clay 

1460 

1.421 

0.1 

0.03  ~  0.2 

Basement 

1800 

2.0 

0.5 

TABLE  I 

Seabed  properties  oe  three  dieferent  sediment  models  and  oe  the  basement.  For  each  sediment  model:  the 

FIRST  LINE  GIVES  THE  OPTIMUM  VALUES  WE  OBTAINED,  THE  SECOND  LINE  INDICATES  THE  DEVIATION  RANGES  OF  THE 
PARAMETER  VALUES  THAT  WE  TESTED  ON.  OVERALL,  MORE  THAN  THIRTY  SEDIMENT  MODEL  VALUES  WERE  COMPARED. 
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Shallow  water 
Shelf  Break 
Continental  Slope 


Z  <  140  m 
140  <  Z  <  210  m 
Z  >  210  m 


Muddy-sand 

Clay 


TABLE  II 

Depth-dependent  Sediment  Model,  resulting  oe  a  best-fit  of  our  coupled  ocean-acoustic  models  to  TL 

OBSERVATIONS.  NOTE  THAT  MORE  THAN  10  SETS  OE  TWO  TRANSITION  DEPTHS  WERE  TESTED,  ONLY  THE  BEST-FIT  IS 

REPORTED  HERE. 
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V.  Figures 


Fig.  1.  a)  QPE  Pilot  Experiment  -  CTD  stations  and  Main  Study  Area.  The  intensive  acoustic  region  is  identified  by  the  solid 
box.  The  solid  dots  are  CTD  stations  for  the  ORl  Leg  1,2  cruises;  stars  for  the  OR2  cruise,  and  circles  for  the  OR3  cruise,  b) 
Bathymetry  (m)  of  the  Acoustic  Experiment  region  overlaid  with  the  two  planned  tracks  of  the  OASIS  Mobile  Acoustic  Source 
(OMAS)  circle  runs,  refereed  to  as  Events  A  and  B.  c)  Bathymetry  (m)  at  locations  of  Events  A  and  B. 
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0"  0" 


Fig.  2.  a)  Measured  mean  TL  of  Events  A  (blue,  with  observed  bearing-averaged  cr  =  1.8  dB)  and  B  (red,  with  a  —  2.4  dB) 
[53].  b)  Coupled  ocean(4D)-acoustics(Nx2D)  simulations  of  the  mean  TL  of  Event  A  (blue,  with  bearing-averaged  a  =  1.23 
dB)  at  12:00  UTC,  September  8,  2008,  and  Event  B  (red,  with  a  =  2.33  dB  in  the  bearing  angle  ranges:  >  260  degree  or  < 
110  degree;  a  =  7.87  dB  for  whole  circle)  at  15:00  UTC,  September  8,  2008. 


January  14,  2010 


DRAFT 


29 


Fig.  3.  Nx2D  TL  simulations  of  Events  A  and  B,  illustrating  effects  of  different  sediments  models  (and  keeping  all  other  factors 
constant:  in  particular,  the  sediment  sickness  is  set  to  20  m).  The  four  different  sediments  discussed  here  were:  depth-dependent 
(hybrid)  model  with  variahle  sediments,  sand,  muddy-sand,  and  silty-clay.  With  these  sediment  variations,  one  obtains  a  standard 
deviation  for  Event  A  of  =  9.55  dB  and  for  Event  B  of  =  5.81  dB.  The  ocean  fields  are 

those  of  12:00  pm,  Sep.  8,  2008  (UTC). 
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(a)  (b) 

Fig.  4.  Nx2D  TL  simulations  of  Event  A  ((TTL^^thymetry  (^)  =  1-25  dB)  and  Event  B  ((^tl bathymetry  W  =  1-02  dB)  with 
different  percentages  of  random  perturbation  on  the  lOOm-resolution  bathymetry  data.  The  normally  distributed  random  noise 
was  added  on  the  bathymetry  with  a  variance  proportional  to  the  slope  and  to  the  local  depth.  The  ocean  fields  are  those  of 
12:00  pm,  Sep.  8,  2008  (UTC). 
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Fig.  5.  Predicted  ocean  sound  speed  field  at  30  m  (a  b  c),  60  m  (d  e  f),  and  90  m  (g  h  i)  depth  at  12:00  pm,  Sep.  8,  2008 
(UTC),  which  is  the  day  when  Events  A  and  B  occurred  at  sea.  Each  column  corresponds  to  a  different  initial  condition  for  the 


vertically-averaged  transport  in  the  Taiwan  Strait:  First  column  -  1  Sv  S;  Second  column  -  0  Sv;  Third  column  -  1  Sv  N. 
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Fig.  6.  Nx2D  TL  simulations  of  Event  A  =  0.58  dB)  and  Event  B  =  1-21  dB)  at  12:00  pm, 

Sep.  8,  2008  (UTC)  corresponding  to  different  ocean  field  predictions  (transport  conditions  between  Taiwan  and  mainland  China 
initialized  at  ±1  or  0  Sv,  as  illustrated  on  Eig.  5)  used  as  inputs  for  the  background  sound  speed  field. 
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Fig.  7.  Predicted  sound-speed  variability  maps  at  a  depth  of  30  meter:  (a)-(i)  correspond  to  the  times  from  00:00  Sep  8  to 
00:00  Sep  9,  2008  (UTC)  with  three  hours  interval.  Panel  (a)  is  the  total  sound-speed  on  00:00  Sep  8  while  Panels  (b)-(i)  are 


differences  between  the  subsequent  3-hourly  sound-speed  fields  and  this  00:00  Sep  8  field. 
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(g) 


(h) 


(i) 


Fig.  8.  Predicted  sound-speed  variability  section,  along  the  shelf,  within  the  acoustic  region:  (a)-(i)  correspond  to  the  times 
from  00:00  Sep  8  to  00:00  Sep  9,  2008  (UTC)  with  three  hours  interval.  Panel  (a)  is  the  total  sound-speed  on  00:00  Sep  8  while 
Panels  (b)-(i)  are  differences  between  the  subsequent  fields  and  this  00:00  Sep  8  field. 
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Fig.  9.  As  Fig.  8,  but  for  the  across-shelf  section  within  the  acoustic  region. 
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(g)  (h)  (i) 

Fig.  10.  Nx2D  TL  simulations  of  Event  A  (blue-curves,  OTL^i^^id)  =  2.21  dB)  and  Event  B  (red-curves,  o-TLti^^{d)  =  3.29 
dB,  in  the  bearing  angle  ranges:  >  260  degree  or  <  110  degree;  a  =  7.64  dB  for  whole  circle  at  different  times).  Panels 
(a)-(i)  again  correspond  to  the  times  from  00:00  Sep  8  to  00:00  Sep  9,  2008  (UTC)  with  three  hours  interval.  This  includes  the 
period  during  which  Events  A  and  B  actually  occurred  at  sea.  Panel  (a)  is  the  total  TL  on  00:00  Sep  8  while  Panels  (b)-(i)  are 
differences  between  the  subsequent  TL  fields  and  this  00:00  Sep  8  field.  The  TL  range  for  these  differences  is  from  -10  dB  to 
10  dB:  the  dashed  circle  denotes  the  0  dB. 
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Fig.  11.  Sections  on  Sep.Sth,  2008,  all  times  UTC:  a)  sound  speed  at  12:00;  b)  sound  speed  at  21:00;  c)  difference  of  sound 
speeds  between  a)  and  b);  d)  TL  estimate  at  12:00;  e)  TL  estimate  at  21:00;  e)  TL  estimates  for  a  receiver  at  61  m,  at  12:00 
(red)  and  21:00  (blue). 
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Fig.  12.  Bearing-averaged  TL  estimates  for  Events  A  and  B  every  three  hours  during  Sep.  8th  (UTC). 
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